Accretion in the spin-down regime 
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Abstract. 

Accretion driven millisecond X-ray pulsars can accrete over a wide range of mass flow rates. The 
pulsations persist even at small accretion rates at which these objects would be expected to be in 
the propeller stage. We argue that the inner regions of disks around millisecond X-ray pulsars are 
sufficiently thick that a fraction of the inflowing matter can accrete even in the spin-down regime 
from regions of the disk away from the disk plane. This allows these systems to be bright throughout 
a wide range of mass flow rates. We model the lightcurve of SAX J 1808.4-3658 during an outburst 
and show that the rapid decay stage can be modeled with fractional accretion in the spin-down 
regime. 
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INTRODUCTION 

The discovery of accretion-powered millisecond X-ray pulsars (AMXP) [1, 2] supports 
the recycling scenario [3, 4] according to which millisecond radio pulsars [5] are de- 
scendants of low mass X-ray binaries. 

All of the 13 AMXPs known to date are X-ray transients. During an outburst the 
luminosity of these objects can change two orders of magnitude. According to the 
conventional picture classifying neutron stars in terms of their interaction with the 
surrounding matter, as the mass flow rate M declines throughout the outburst, these 
objects should undergo a transition to the propeller stage [6] in which the centrifugal 
barrier inhibits inflowing matter in the disk to reach the surface of the neutron star. 
Around this transition the torque should also change sign so that the star spins down. 

Accretion onto AMXPs even at low accretion rates is a challenge to this picture. In 
some cases, even while the the neutron stars exhibit a high rate of spin-down, accretion 
onto the neutron star continues, and the X-ray pulsations prevail. In order to explain 
these difficulties Rappaport et al. [7] has presented a steady state solution of the disk 
structure in which the inner radius of the disk remains on the co-rotation radius while M 
declines. 

SAX J1808. 4-3658 experiences outbursts lasting a few weeks roughly once in two 
years, during which the coherent ~ 401 Hz pulsations are observable. The inner radius 
of the accretion disk, R m = t, (/l 2 / y/lGM^M) 2 ! 1 where /! is the magnetic moment, M* 
is the mass of the neutron star and £ is a factor of unity, is expected to increase with 
dropping M. There is recent evidence [8, 9] that the R m does recede during the decline 
of the outburst. 

An infinitely thin disk can not accrete matter as soon as i?i n moves beyond the co- 
rotation radius R c = (GM*/^) 1 / 3 where Q.* is the angular frequency of the neutron 



star. In this case the fastness parameter ft)* = (R- m /Rc) 3 ^ 2 becomes greater than unity. The 
luminosity is expected to drop to zero abruptly. The situation can be different for a disk 
with a. finite thickness which will be able to accrete matter even for ft)* > 1 because the 
part of the disk at high latitudes will not be centrifugally inhibited. Inclination between 
the rotation axis and the magnetic axis will further facilitate accretion. 

Here, we reassess the work of [10] for estimating the fraction of mass that can accrete 
in the spin down regime and apply it to model the lightcurve of SAX J1808.4-3658. 



MATTER ACCRETING IN THE SPIN-DOWN REGIME 

The fraction of matter that can accrete onto a neutron star in the spin-down regime 
was first estimated by Lipunov & Shakura [11] for spherical accretion. Here we follow 
the work of Menou et al. [10] for quasi-spherical disk accretion [14]. There are recent 
numerical studies [12, 13] that address accretion in the spin-down regime. 

The magnetic surface is defined as r m = Rm(l + 3cos 2 0) 2//? where measures the 
angle from the rotation axis which we assume to coincide with the magnetic axis. The 
co-rotation surface is defined as r c = R c sin -2 / 3 6. A disk with a finite thickness will be 
able to accrete matter even when Ri n > R c because at high latitudes (0 < < Oo) r m < r c 
can be satisfied. The fraction of the inflowing matter that can accrete onto the neutron 
star is 

_M*_ 2 f°° pvlnr 1 sine d6 

~ M 2f* /2 pv2Kr 2 sm9de' 

Here the critical angle 0q is defined through r^ = r c . Given the uncertainty of the size of 
the magnetosphere in the presence of plasma, we follow [10] and assume that r m ~ R m 
which gives sinOo = ft)" 1 . Using p(r, 6) ~ p(r) and v(r, 6) ~ v(r) sin 2 from the quasi- 
spherical solutions of [14], one obtains 

= 1 - 1 (1 - co- 2 ) 1/2 + 1 (1 - co- 2 ) 3/2 (2) 

Note that this result follows from the work of [10], the only difference being that these 
authors, at the final step, used the small angle approximation (i.e. 0q<^\) and so found 
/ = (3/8) ft)" 4 which is ft)* > 1 limit of the above. 



EVOLUTION OF THE DISK 

The thick inner region of the disk will match to a standard thin disk [15] outside. The 
evolution of M during the outburst will be determined by what the outer thin disk 
provides to the inner quasi-spherical region. The evolution of a thin disk is described 
by the diffusion equation 
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FIGURE 1. The lightcurve of SAX J1808.4-3658 in October-November 2002 outburst compared with 
the model presented in this work. The data is taken from the first figure in [8] excluding the flaring decay 
stage which is not addressed here. The peak stage can not be modeled with the self-similar model for M 
employed in this work. The slow decay stage is modeled by assuming all the inflowing mass accretes onto 
the star. As M declines o* grows beyond unity. In the rapid decay stage only a fraction of M can accrete. 



where E is the surface mass density and v is the turbulent viscosity. For a standard thin 
disk [15], this equation has three self-similar solutions [16] two of which have freely 
expanding outer radius and are not suitable for tidally truncated disks in binaries. The 
third solution is suitable for disks in binary systems [17]. According to this solution 



where Mq is the mass flow rate at the beginning of the outburst, to is the viscous time- 
scale in the outer disk, and a depends on the opacity and pressure regime prevailing in 
the disk. For an electron scattering and gas pressure dominated disk a = 5/2. 

The X-ray luminosity of accretion onto the neutron star is found by Lx = GMM*/R 
where M* = M if co* < 1 and M* = fM if to* > 1 . 



SAX Jl 808.4-3658 showed an outburst in October-November 2002 followed by Rossi 
X-ray Timing Explorer. The lightcurve shows 4 stages: peak, slow decay, rapid decay, 
and flaring decay as described by Ibragimov & Poutanen [8]. We converted the 3 — 20 
keV flux data therein to luminosity by assuming the distance to be 3.5 kpc [18]. 

The slow decay region can be modeled by using to = 19 days, Mo = 2.2 x 10 16 g s _1 
and a = 5/2 in Eqn.(4) assuming all matter can acccrete. 




(4) 



APPLICATION TO SAX J1808.4-3658 



We associate the commence of the rapid decay with ft)* = 1 which allows for deter- 
mining the magnetic moment of the star as jl = 0.7^ ~ 7 / 4 x 10 26 G cm 3 . For £ = 0.5 
this corresponds to a magnetic moment of /i = 2.2 x 10 26 G cm 3 . We assume that M 
declines at the same rate given in Eqn. (4), but now a fraction, determined by Eqn.(2), 
of it accretes onto the neutron star. We find that the rapid decay stage can be modeled as 
fractional accretion in the fast (ft)* > 1) regime. 

Because of the self-similar model employed for M, the rise to the peak stage cannot 
be addressed in this work. This can be done by solving Eqn. (3) numerically with an 
appropriate initial distribution for E. The rapid decay stage is followed by a flaring decay 
stage (see the Figure 1 in [8]) which is not addressed in this work. 

DISCUSSION 

We have suggested that the rapid decay stage of the lightcurve of SAX J 1808.4-3658 
during an outburst can be modeled by fractional accretion in the rapidly rotating stage 
(ft)* > 1). This requires the inner region of the disk to be sufficiently thick. Our result, 
B = 1.4^ -7 / 4 x 10 8 G, is consistent with earlier estimates [19, 20, 21]. 
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